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Abstract 

Starting with a percolation model in Z d in the subcritical regime, 
we consider a random walk described as follows: the probability of 
transition from x to y is proportional to some function / of the size 
of the cluster of y. This function is supposed to be increasing, so that 
the random walk is attracted by bigger clusters. For f(t) = e^* we 
prove that there is a phase transition in (3, i.e., the random walk is 
sub diffusive for large (5 and is diffusive for small (5. 

Keywords: subcritical percolation, subdiffusivity, reversibility, spec- 
tral gap 



1 Introduction and results 

First, we describe the usual site percolation model in Z d . This model is 
denned as follows. For fixed p G (0, 1), consider i.i.d. random variables u(x), 
x G Z d , where u(x) = 1 with probability p and u(x) = with probability 
1 — p. A site x is said to be open if u(x) = 1 and closed otherwise. Write 
x ~ z if x and z are neighbors. A (self-avoiding) path from x to y is: 
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j(x, y) = {xq = x, Xi, x 2 , ■ ■ ■ , x n = y}, where Xi ^ Xj if i ^ j and X; L ~ x i+ \, 
% = 0, . . . , n — 1. A path 7 is said to be open if all the sites in 7 are open. 
The cluster of x is defined by 

<£(x) = {y £ Z, d : uo(y) = 1 and there is an open path 7(0;, y) from x to y}. 

Note that, if oj(x) = 0, then €(x) = 0. It is a well-known fact (see e.g. jU]) 
that there exists p cr (depending on d; obviously, p cr = 1 in dimension 1) such 
that if p < p cr , then a.s. there is no infinite open cluster, and if p > p cr , 
then a.s. there exists an infinite open cluster (also, with positive probability 
|C(0)| = oo). 

Throughout this paper we assume that p < p cr , i.e, the model is in 
the (strictly) subcritical regime. Fix a parameter j3 > 0. The percolation 
configuration is regarded as random environment. Fixed the environment, 
we start a discrete time random walk on Z d with transition probabilities 

e mv)\ 

pu) _ 

if x ~ y. Since j3 is positive, one can note that the random walk is in some 
sense "attracted" by bigger clusters, and the strength of this attraction grows 
with (3. Denote by £ (t) the position of this random walk at time t. Let P be 
the probability measure with respect to to and the (so-called quenched) 
probability for the random walk starting from x in the fixed environment uj. 
Denote also P x = <S> P (P x is usually called the annealed probability); 
throughout the paper || • || stands for the norm. Our main result is 
that there is a phase transition in j3, i.e., the random walk exhibits different 
behaviors for large and small (3: it is diffusive for small values of (3 and 
subdiffusive for large values of (3. 

Theorem 1.1 Suppose that the random walk £(t) starts from the origin. 
There exist /3 and (3\ (depending on d) such that < /3 < (3\ < 00 and 

(i) if0</3< (3 , then 

lim logmaxo< s < t il^)H = 1 p0 _ a s 
t^oo logt 2 



(ii) if P > Pi, then 



logmax 0<s <t ||£(s)|| 1 

hmsup 7=-^ < -, P -a.s. 1.2 

t^oo logt 2' v ; 
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One can prove also that the same result holds for the bond percolation 
model in the subcritical regime. The method of the proof remains the same; 
the reason why we have chosen the site percolation is that for bond perco- 
lation there are some technical difficulties (easily manageable, though; they 
relate to the fact that, in the bond percolation model, two neighboring sites 
can belong to different large clusters) in the proof of the part (ii) of Theo- 
rem 

Recently much work has been done on the (simple or not) random walk on 
the unique infinite cluster for the supercritical (bond or site) percolation in 
Z d (see e.g. j3 IH HH1 US] ; see also [TT] for some results for the random walk on 
the incipient infinite cluster in dimension 2). Another related subject is the 
class of models (see e.g. 0E]) that can be described as follows. Into each edge 
of Z d we place a random variable that represents the transition rate between 
the sites. The new features of the model of the present paper are, first, 
the fact that the random environment is not independent, and secondly, the 
absence of the uniform ellipticity. Speaking of uniform ellipticity, we should 
mention that in the paper [7j there was considered a simple symmetric one- 
dimensional random walk with random rates, where the time spent at site i 
before taking a step has an exponential distribution with mean 7$, and Tj's 
are i.i.d. positive random variables with distribution function F having a 
polinomial tail. One may find that there are similarities of the <i-dimensional 
analog of the model of |X with our model, because clusters of size n will have 
"density" e~ Cn , and the mean time spent there is roughly e 1371 , so, thinking 
of clusters as "sites" , we indeed obtain a polinomial tail of mean time spent 
at a given site. However, the facts that the random environment is no longer 
independent and that here the random walk is not a time-change of the 
simple random walk make the model of the present paper considerably more 
difficult to analyze. 

2 Proof of Theorem 11.11 

We begin by introducing some notations and recalling a well-known fact from 
the percolation theory. Namely, we will use the following result (see |12| 19]): 
if p < p cr , then there exists c\ > such that for all N > and 

P[|£(x)| > N] < e~ ClN . (2.1) 
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Now, to prove Theorem 11.11 an important idea is to consider £ (t) in finite 
region. Take A n = (— n/2,n/2] d and let the process ^ n '(t) be the random 
walk restricted on A n . 

Proof of part (i). It can be easily seen that £ (£) is reversible with the reversible 
measure 

n(x) = e^Y, emZ)l i ( 2 - 2 ) 



and thus the finite Markov chain ^ n '(t) is also reversible, with the invariant 
(and reversible) measure 

e p\<t{*)\ e p ^\ 
*<»>(*) = z -f , (2.3) 

where 

is the normalizing constant, so that XLeA ^ i x ) = 1- 

Consider also a random walk £^(t) that is a continuization of ^ n \t). 
That is, = £ {n) (iV t ), where iV t is a Poisson process with rate 1, indepen- 

dent of anything else (in other words, is a continuous time Markov chain 
with the transition rates equal to the transition probabilities of £ <n - ) ). Let % 
be the time interval between the jumps (i — 1) and i of N t , S n = X^=i^' 
and Ta (respectively, Ta) be hitting time of set A by random walk ^ n \t) 
(respectively, C (t))- It can be easily seen (cf., for example, Chapter 2 of [T]) 
that Ta = St a and E(T^ | Ta) = Ta- Moreover, since t -1 iVt — > 1 a.s., many 
other results concerning £( n ) can be easily translated into the corresponding 
results for 

Remark 2.1 Using this technique, it is elementary to obtain that Theo- 
rem holds for £(£) iff it holds for £(t), where £(t) is the continuization 
°f £{t) defined in the same way. 

So, now we consider the finite continuous time Markov chain ^ n \t). De- 
note by A the spectral gap of £ ft)- 

Lemma 2.1 There exist C14 > and n* = n*{u) such that for all n > n* we 
have A > c u n~ 2 , F-a.s. 
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Proof of Lemma \2. 1\ The idea is to use Theorem 3.2.1 from (Hj to prove 
the lemma. For each pair x, y G A n , we will choose exactly one path 7(x, y) 
(connecting x and y) in a way described below. Let \^{x,y) \ be the length 
of 7(2,?/) (i.e. the number of edges in 7(2,;?/)). Denote by £(A n ) the set 
of edges of Z d n A„. For an edge u = (z lj z 2 ) let Q{u) = (P^^ (zi) + 
P^ 2z y n \z 2 ))/2. According to Theorem 3.2.1 of [H], it holds that A > l/A, 
where 

A= max {J- £ |7(a:,y)k (B) (^ (B) (y)}- ( 2 - 4 ) 
«e£(A w ) I Q(u) „, , > 

x,y£A n : 7(x,j/)9m 

Here, we have, for u = (zi, z 2 ), 

e^(|«(*i)|+|(£(«a)|) 
Q(u) = = 



and for each pair x = (x^\ . . . , x^),y = (y^ 1 *, . . . , y^) G A n we choose the 
path j(x,y) in the following way. Let ei,...,e d be the coordinate vectors. 
Denote A, = g/W — x^' and let sgn(Aj) be the sign of Aj. Suppose for 
definiteness that x^ < j/W (so that sgn(Ad) > 0). We take then 

~/(x,y) = (x,x + e d , . . . ,x + A d e d ,x + sgn(A d _i)e d _i + A d e d , . . . , 
x + A d _ x e d ^ x + A d e d , ...,x + Aid H 1- A d e d = y), 

so, first we successively change the <i-th coordinate of x to obtain the d-th 
coordinate of y, then we do the same with (d — l)-th coordinate, and so on. 
With this construction it is clear that the length of 7(2, y) is at most dn. For 
an edge 

u= ((xW 1 ...,iW),(x (1) ) ... ) x (d - 1) ,x (d) + l)) 

define 

4 = {(*« . . . , *«>) G A„ : = jd), . . . , = a M < S M} 



and 



i? u = {(z (1) ,...,z {d) ) G A n : 2 W>a;W} 



(for the edges of other directions the computations are quite analogous). We 
have then 



E 



7(x,y)KW(x)7rW(y) 
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xel u y&Ru 



(2.5) 



Now our goal is to prove that with large probability, for all such u, 

X&I U 

Iu = {{z {1) ., 



^2 x€l 7c^ n \x) is of order n/Z. Denote 



z W)eK: z<n=xM,...,z( d -v=x«- iy >}. 



Note that I u C I u , so we will concentrate on obtaining an upper bound for 
^ xg j 7[( n \x). It is important to observe that the variables ix^ n \x) are not 
independent. For the sake of simplicity, suppose that y/n is an integer, the 
general case can be treated analogously. Divide I u into sfn equal (connected) 
parts of size y/n, denote := (x^\ . . . , x^ d ~ l \ (i — l)\/n + j) and write 

I — Vn 
y/n — 

E- (n) (*) = E E ^W- 



-+1 



For fixed j, let 



L{| I <>/S/2}> 



and B = DiS- 8 *- % dB> we have 
> 1 — e~ Clv/ ™ which implies that P[5] > 1 — e~ C2v/ ™ for some C2 > 0, so 

k\/n~ 



2 



P 



[ £ 



< P 



y n 
2 



E ^ (n) o%)> 



lb 



+ e 



y/n 



+1 



Now, it is important to note that the variables e^' £ ^ lJ '-" Y^ z ~ x - e^'^'ls = 
7r(xy)ls (recall (|2.2|) and (|2.3|l ). z = 1, . . . , \/n, are independent. We have 
also 



E < 2rfe 



2 J S|£(x ij -)l 



where x^- satisfies |<£(%)| = max^^C^-)!, |<£(z)|}. 
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For y = n a (a > will be chosen later), using (|2.1j) . we have 

W[2de 2m ^ )l >y]< v / ™r £a ^ = v^?/~i = n*~^. (2.6) 

According to Corollary 1.5 from [13] . if Xi, . . . ,X k are independent ran- 
dom variables, Sk = Yli=i Xi, Fi{x) = ¥[Xi < x], then for any set y\, . . . , y^ 
of positive numbers and any t, t e (0, 1], 

eAt 



nsk >x]< J2nx > Vi] + (^i) ? > ( 2 - 7 ) 

where y > maxj?/!, . . . , y k } and 

k 00 

At = J2 [ u'dFiiu). 
1=1 o 

Denote Fi(u) = 1 — Fiiu). We apply Corollary 1.5 from ^3] to random 
variables 2de 2 ^ € ^l B , i = + 1, . . . , ^, with x = ky/E, y { = y = n a , 
and t — 1. First term of the right-hand side of ([2.7]) was estimated in ([2.6]) . 
For the second term, we have 



\A> oo 
2 



At = J2 j udF * 



u 



i=-^+l 



OO 

— \fn J udF L {u) 
o 

oo 

— \pa J Fi(u)du 



o 

oo 



< c^yfn I u ^du 



1 



csv 7 ^, (2-8) 
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as (3 is small, thus 



/ — fcv/n 

y < i ec 5 y/n\ — 



k^fn 



(t) 



i 



(2.9) 



so, to guarantee that 



eA+ 



as n — > oo, it is sufficient to take k large 



enough and /3/c 3 < a < 1/2. 

We proved that for (3 sufficiently small 



2 



P 



< c 7 n 2{ p ' < c 7 n + 2 . 



(2.10) 



Thus, 



P [5Z 7r(n) (a ; ) > 



Csn 
Z 



V™ — 



P 



[E E - (n) (^)> 



< 



E p E 



i=i 



c 8 n 



ky/n 



< c 7 n + 



So, using (|2~TH) in (|2~Hj) and (|2~l|l . we have 



(2.11) 



A < max 



c 9 Z 



rr 



ue£(A„) l e /3(l^i)l+l^2)|) ^ 



< CgH 



[2.12) 



as e^d^^^l+l^ 22 )!) > 1, with probability at least 1 - c 7 n -f ^ +1 . Since /3 can 
be made arbitrarily small, Borel-Cantelli lemma implies that for almost all 
environments for n large enough it holds that A < c^n 2 and thus A > C\^rT 2 . 
Lemma f2. II is proved. ■ 

Now, using Lemma 2.1.4 from [TJ] with f{x) = l{|| a; ||>n/4} ) where, as 
before, || • || is the norm, we prove (jl.lj) . By Lemma 2.1.4 from ^1] we 
have that 

\\H t f-n( n \f)\\l<e- 2Xt V^\f). 
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In what follows we show that ^ n \f) is of constant order. We have 

e /3|£(x)| ^ e^ 2 )! 

^(/) = E — z -f — ■ 

n, ||x||>7l/4 

Since [€(a;)| > 0, it is easy to obtain that for all u; it holds 

J2 e m*)\J2e^>-. (2.13) 

Using the same kind of argument as in the proof of Lemma I2~T1 one can easily 
see that for all n 



P 



e *Wl^ e »)l> Cl / <c'>-^\ (2.14) 

xeA n , \\x\\<n/4 z~x 



where c' 15 , c" 5 depend only on C15. Thus, with probability at least 1 — c'{ 5 n~ p 
we have n^ n \f) > const. Then, using that Var^(/) < 1, taking t = ciqu 2 
for Ci6 large enough yields that the random walk ^ n \t), and thus ^ n \t), will 
be at distance of order n from the origin (as both random walks start from 
0) after a time of order n 2 with probability bounded away from 0. 

Now, for any fixed e > 0, divide the time interval (0, t] into t e intervals of 
length t 1_e . Borel-Cantelli lemma implies then that for t large enough there 
will be at least one time interval such that at the end of this interval 
will be at distance at least t^~2 from the origin. Since e > is arbitrary, we 
proved that 

liminf logmax <^||^)|| > 1 o 
t^oo log* ~ 2 v 1 

It remains to prove that 

logmaxo< s<i ||£(s)|| 1 

hmsup < -, P -a.s. (2-16) 

t^oo log t 2 

It is a well-known fact that a reversible Markov chain with a "well-behaved" 
reversible measure cannot go much farther than t 1 ^ 2 by time t, see |3J El EU 
ITd] . By Theorem 1 from jB] we have, for any e > 

2(l/2+ e ) P PWV)\ V p/W)l 

^liiuii > n j s ^ 2^ e mm V e^ 2 )! 



9 



< c 20 e- nE P-a.s. 



(2.17) 



for all n large enough. To obtain the bound (|2.17|) we have used the fact 
that, due to (EHj) . 

F\maxe mx)l > n*} < n'^ 

for some C21 > 0. Borel-Cantelli lemma and ()2.17|) imply (|2.16|) and thus the 
part (i) of Theorem 11.11 is proved. 

Proof of part (ii). For iGZ d let T (x) = 0, Tq(x) = 0, and define 

T'^x) = min{i > T^x) + Tt_ x {x) : £(t) G £(x)} 
^(rr) = min{t > ^'(x) : f (*) g <£(a:)} - Tftz), 

i = 1,2,3,..., where Tj(x) is defined if min{t > T k (x) : £(£) G £(x)}, 
k < i, are finite. In words, T/(x) is the moment of ith entry to the cluster 
of x, and Tj(x) is the time spent there (i.e., after T/(x) and before going 
out of C(x)). It is important to note that the cluster C(x) is surrounded by 
sites with uj(-) = 0. Comparing T,(x) with geometric random variable with 
parameter (2d — 1 +e^ c I 0')~ :L , one can easily see that if > 51ogn, then 

ETj(x) > cirra^ 5 . Moreover, it is elementary to obtain that for e > and for 
any 5 > we can choose /3 large enough so that with probability bounded 
away from we have 

Ti(x) > c 18 n 2+£ (2.18) 

for x such that \C(x)\ > 5\ogn. 

Now, we use a dynamic construction of the percolation environment usu- 
ally called the generation method (see [12] )• That is, we proceed in the 
following way: we assign generation index to the origin, and put 00 (0) = 1 
(the origin is open) with probability p or a>(0) = (closed) with probability 
1 — p. If co>(0) = 0, the process stops. If lu(0) = 1, then to all x ~ we 
assign generation index 1, and put u>(x) = 1 with probability p or w(x) = 
with probability 1 — p, independently. Suppose that the we constructed m 
generations of the process. Let Y; be the set of sites with generation index % 
and Y m = {0} U Y x U . . . U Y m . Denote by Y m+1 the set of neighbors of the 
open sites in Y m which do not belong to Y m . Assign to the sites from Y m+ i 
the generation index m + 1 and a value 1 or in a way described above. If 
Y m 7^ and u)(y) = for all y G Y m+ x, then the process stops. Note that 
for subcritical percolation this process stops a.s., and what we obtain at the 
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moment when the process stops is the cluster of the origin surrounded by 
0-s. 

So, first we construct the environment within the set H\ = £(0) U <9£(0), 
where 

<9£(0) = {y : y (£(0), y ~ x for some a; G £(0)} 

(note that uj(y) = for any y G <9£(0)) and we know nothing yet about the 
environment out of the set H\. For an arbitrary set H C Z d denote 

H° = {x H : for any infinite path 7(2) starting from x it holds 
that 7(ac) n H ^ 0} 

(i.e., if° is the set of the holes within the set H) and let 

G l = H 1 U H°. 

Then, choose u(x) for x G H° and start the random walk £(t) from the origin. 
Let 

n = min{t : f (t) £ G?i}. 
Note that £(£(ti)) fl £(0) = 0, and construct, using the above method 

F 2 = C(e(r 1 ))u9C(e(r 1 )) 

and 

G 2 = d U # 2 U (G 1 U tf 2 )°. 

Then, define 

r 2 = min{t : £(t) t G 2 }, 
and so on. For all i, we have 

P[C(f(Ti)) > 51ogn] >p 51 °s" 

where p is the percolation parameter. This is so due to the fact that, to have 
C(£(Ti)) > 51ogn, it is sufficient to choose a path of length <51ogn emanating 
from £(11) which does not intersect Gi (it is possible by the construction of 
Ti, since £(Tj) cannot be completely surrounded by points of Gi), and such 
path will be open with probability p si °g n . Thus, for any e > (one can take 
the same e from (|2.18|) ). we can choose 5 small enough (take 5 such that 
<51ogp -1 < e) so that 

P°[£(^)) > ^ogn I Ti] > n- £ , (2.19) 
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where Ti is the a-algebra generated by {u)(x), x G Gi} and {£{m), m < Tj}. 
Fix 8 > in such a way that 1 — 9 > e. Note that, as p < p cr , using (|2.1|) 
and Borel-Cantelli lemma, for n large enough min{fc : £(Tfc) ^ A n } (the 
number of times that we repeat the basic step in the above construction) 
will be of order at least n 1 ^ 9 for all n large enough, P°-a.s. (recall that 
A n = (— n/2, n/2] d ; with overwhelming probability all the clusters inside 
A n will be of sizes at most n e ). On each step, by (J2.19|) . with probability at 
least n~ £ the random walk enters the cluster of size at least 5 log n. By ()2.18|) . 
it stays in that cluster (if f3 is large enough) for at least n 2+e time units with 
large probability. If 1 — 9 > e, with overwhelming probability on some step 
(of the above construction) the random walk will delay (in the corresponding 
cluster) for more than n 2+e time units before going out of A n . In other words, 
we will have 

max ||£( s )||<^, 

which implies (11.2)1 . This concludes the proof of Theorem 11.11 ■ 
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